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Abstract 


Great effort is being focused on making the next generation of naval combatant ships 
more resistant to the effects of close-aboard explosions. The examination of the 
deformation modes in blast-loaded metal plating suggests that a physical model can be 
developed to simulate the force vs. displacement history produced by an impinging shock 
wave during the holing phase. Similar approaches have been successfully used to 
approximate damage due to grounding and ballistic penetrators. 

In this case, the deformation of the clamped plate is modeled in two stages: (1) dishing, 
which leads to disking and (2) radial crack propagation, which results in petalling. In the first 
stage, a thin geometrically-scaled (0.90 mm, 1.15 mm, and 1.40 mm thick by 300 mm square) 
mild steel sheets are dished inward using spherical indenters of radii 20 mm, 50mm, and 75 
mm. The sheets have an average tensile strength of 317 MPa and a Rockwell Superficial 
Hardness Number of 72 (H N15T 72). 

This portion of the test approximates the initial material stretching done by a spherical 
wave at various standoff distances. The spherical indenter produces a circular hole, which 
simulates the disk of material normally ejected as a blast front penetrates a plate section. As 
the material reaches a critical necking thickness at the edges of the hole, radial cracks form 
creating petals. During the second stage, an oblique conical punch is used to simulate the 
expanding wave front, which drives open the petals, causing the cracks to propagate towards 
the plate’s clamped boundaries. By measuring the resultant forces and minimizing the 
effects of friction, the total bending and membrane work can be reasonably estimated. 
Ultimately, the approximate blast damage for a given ship’s hull may be related to a given 
charge size and standoff distance. 
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Introduction 


In response to the growing number of threat nations and the increased proliferation of 
anti-ship weaponry, the Office of Naval Research (ONR) has prompted new research in the 
areas of weapons effects and ship vulnerability. Mindful of budgetary pressures, the United 
States Navy is exploring the most cost-effective methods of increasing the damage resistance 
and improving the overall battle effectiveness of its warships. (Refer to Appendix A for 
more detailed prefatory material.) 

To date, the plastic deformation resulting from a close-proximity explosion (either above 
or below the waterline) has not been the focus of any large-scale, publicly-accessible research 
project. Many organizations have intensely researched topics which apply similar mechanics, 
such as ballistic penetration, collision and grounding damage, and axial tube splitting. There 
has been preliminary work by both Wierzbicki (1996)/(1999) and Nurick (1996), which 
examined plate tearing and petalling done by on contact explosives. (A detailed review of 
past work in this area and the associated literature can be found in Appendix B) 

In keeping with the goals of the U.S. Navy and ONR, this research focuses on the 
deformation and fracture of hull plating, which is subjected to either an underwater 
(UNDEX) or air explosion. Ultimately, the objective of the study is to provide a simplified 
method of gathering benchmark data to quantify a material’s sensitivity to explosive damage. 
An inexpensive two-stage quasi-static indentation test is used to model the force vs. 
displacement history of a thin clamped steel sheet. These experimental results are then 
compared to both approximate analytical solutions and results obtained through Numerical 
simulations created in ABAQUS (Static Load with Pressure Boundary Conditions) and LS- 
DYNA (Quasi-static Load with Displacement Boundary Conditions). 
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Nomenclature 


R Plate Width 

r Hole Radius 

Rt, Spherical Punch Radius 

Rc Clamped Plate Width 

F f Frictional Force 

P Indenter Force 

V Indenter Velocity 

p Frictional Coefficient 

o y Material Yield Strength 

% Material Ultimate Strength 

a Q Material Flow Stress 

8 Plate Central Displacement 

6 Cone Punch Angle 

<p c Indenter Wrapping Angle 

t Plate Thickness 

E Indenter Work 

1 Petal Length 




Statement of Problem 


To improve blast damage resistance in future classes of combatant ships, the U. S. Navy 
and ONR are entertaining new concepts in both design and materials. Although, the 
characteristics of traditional mild and high strength steels have been studied extensively in 
naval applications, no simple, reliable method exists to predict hull panel blast damage due to 
a close proximity blast load. Blast holing prediction is a key factor in assessing ship 
survivability, including the number of flooded compartments, the ships residual section 
modulus, and probability of recovery. Blast damage computer codes are plentiful, but the 
accuracy of their results often remains in question. Consequently, nearly all trusted damage 
prediction is done through scaled live fire testing, which is both expensive and time 
consuming. 


DISHED SECTION 



Figure 1: Simplified Blast Loaded Plate Geometry 
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The long-term challenge of this research is to combine theories from tearing fracture, 
plate cutting, petalling, and blast loading so as to provide a preliminary estimate of the 
damage caused by a given charge at some small distance from a clamped metal plate. 

The simplified physical model in Figure 1, presents the geometiy of the problem. Figure 2 
shows clear deformation and fracture similarities between dynamic (Rajendran (2000)) and 
quasi-static test specimens. 



On Contact Blast Test (20 g PEK-1 explosive) Quasi-statically Indented Plate 


Figure 2: Dynamically and Quasi-statically Fractured Specimens 

The quasi-static tests have two immediate purposes. First, if the time-pressure history 
of an explosive can be related to the force displacement history .of the proposed quasi-static 
tests, it becomes possible to relate the energy of the incident blast wave to the bending, 
stretching, and tearing work done in the material. Secondly, by examining crack initiation 
and propagation and measuring incremental strains during these quasi-static tests, one can 
develop a rudimentary crack-propagation criterion. Such a criterion could then be used to 
improve existing finite element codes and coupled hydro-codes. It is believed that these 
tools will lead to improve damage estimates (e.g. resultant hole size for a given charge). If 




such estimates can be made quickly and inexpensively, research facilities could rapidly 
examine the suitability of a wide range of materials and structural arrangements with minimal 




Formulation of Problem 

Experimental Approach 

A series of penetrator tests is used to model the deformation which occurs in mild steel 
panels due to close proximity explosions. The modeling occurs in two stages. During the 
first stage, forces, displacements, and local strains are measured as a hemispherical punch is 
used to dish the center of a clamped thin mild steel sheet. As the material stretches and 
thins, a critical necking thickness will be reached near the center of the sheet. A 
circumferential crack then forms at the necked location. Radial cracks will subsequently form 

at the newly formed hole’s edge. The initiation of radial cracks will mark the conclusion of 
the first stage. 

In the second stage, an oblique conical punch is used to examine the behavior of the 
material while petalling. The cone is used to propagate either pre-cut or naturally formed 
radial cracks, causing several distinct outward-opening petals to form. With Teflon-based 
dry film lubricants baked onto both the test specimens and indenter contact surfaces, 
friction is kept to a minimum during the test. Fine grid markings (2.5 mm squares) are 
made on the plate in order to compute final local strains at periodic locations along the 
plate s radius. These strains can than be used to measure material stretching, petal 
displacement/curvature, bending work for the final deformed geometry. This total work can 
than be loosely related to the energy released during the first several milliseconds of an 
explosion. 

Additionally, a series of successive photographs captures the formation and propagation 
of radial cracks through the mesh. In doing so, one can estimate the incremental strain on a 
given mesh element and estimate Crack Tip Opening Angle (CTOA) throughout the 
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propagation phase. When joined with a known force-displacement history, one can begin to 
draw conclusions about the material’s ability to resist blast damage. 

Numerical Approach 

In order to gain insight into the fracture and displacements produced expected the quasi¬ 
static dishing experiment, a numerical model was created in ABAQUS. An axi-symmetric 
plate model was statically loaded with full-plate and half-radius pressures of varying 
magnitude. This analysis showed the (i) mode and location of fracture, (ii) whether fracture 
occurs buy thinning or necking, and (iii) an approximate deformed plate shape at fracture. 


q * 

$ T.~- 

OH] IT 1 ITTT1_5 

3.t 


a t b 
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Figure 3: ABAQUS Numerical Model Representation 

Additionally, a moire comprehensive (and considerably more accurate) model was 
developed in LS-DYNA using displacement boundary conditions for a spherical punch. 

Approximate Analytical Solutions 

Approximate analytical methods have been previously developed for each of the quasi¬ 
static phases. Simonsen (2000) offers an approximate closed-form solution (under several 
assumptions) for plate displacement under hemispherical punch loading. With that solution, 
strains and dishing work can be computed with reasonable accuracy. These results are then 
compared with those obtain experimentally and numerically. 
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Figure 4: Generalized Force vs. Displacement History 

With a closed-form solution derived originally for plate cutting, Wierzbicki (1993)/(1999) 
developed expressions for bending and membrane work done during the circumferential 
cracking/petalling phase. Combined, the various methods provide a reasonable estimate of 
the total work done in deforming the plate during a close-proximity explosion. 

Still, knowing the work of macro-scale deformation is less important than identifying the 
overall force-displacement history. The deformed plate is the most obvious physical clue in 

i 

surmising the wave front-plate interaction from the time of impart through to the petalling 
phase. It is this load-deformation relationship, combined with a suitable crack propagation 
criterion, which can most closely link the quasi-static conditions of a punch test to the 
dynamic conditions present in an explosive blast. 







Approximate Analysis 

Hemispherical Dishing and Circumferential Cracking Phase 

As mentioned earlier, the dishing analysis is based on work done by Simonsen. The 
theory assumes an axi-symmetric plate is dished by a hemispherical punch of radius R b . A 
Cartesian (w, r) coordinate system is used with tp c being the angle from the center of the 
punch to the outermost contact point C. The total punch displacement is 6 and the punch 
force is P. Figure 5 illustrates the geometry of the problem. 



Figure 5: Spherical Dishing Geometry 

Several assumptions are made in order to reach a closed-form solution. Plate bending is 
neglected, plate elements are considered to be displaced vertically, and the material is 
assumed to be rigid-plastic. Under the assumption of plane strain, the generalized radial 
membrane force is given by 





The punch force (P) at point C can be expressed as: 

P = 2 7tN 0 R b sin 2 y/(r) 

where y/(r) is the wrapping angle function 


(2) 


Given the relation that: 


dw 

sm y/ = — 

Jdw 2 +dr 2 


( 3 ) 


The following differential equation can be surmised: 

dw _ sin 2 ^ f 
dr 1/ V 


1 


r 


-sin 4 y/ c 


( 4 ) 


By applying the boundary condition w{r) = 0 at r = R, a solution for w(r) can be reached 
by separation of variables: 


Mr) = 


p 

2/tNo 



r + -Jr^-shT 4 ^" 
R + >//? 2 -sin 4 y^ c 


for r c < r < R 


( 5 ) 


Consequently, using the previously stated assumptions, one can determine the vertical 
displacement for all points on the plate and the force required to deform the plate. Of 
course, at some point the plate the plate will reach its material limits. Strain hardening 
characteristics must be introduced and the true stress strain curve is assumed to obey a 
power law: 



To predict necking localization and fracture, the maximum on the load-displacement 


dp 

curve is identified by setting- 

dy/ c 


= 0 . One can then solve for an approximate solution for 


<p c at necking failure. The final expression given by Simonsen is: 


¥c_jm *-957 + .399h (7) 

Using this relation, additional expressions were derived to determine total displacement (S fai] ) 
and total work (E) up to the point of necking failure: 

6 M ,=\Aln*R M R* (8) 


.318 — 



. 607 - 387 —+ 1.2 
R 


(£)' 


E = nC 0 RR b \ 


+ .067(n-.2) 


( 9 ) 




Conical Dishing 

In some cases, a cone was used to dish a flat plate with a pre-cut starter hole in order to 
gain insight as when and to how the radial cracks form and propagate. 



Figure 6: Conical Dishing Experiment 
The problem’s geometry is defined in Figure 7. 



Figure 7: Conical Dishing Geometry 
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Again, the material is assumed to obey a true stress power law. 
a = CJ" 


(10) 


Using the theory of moderately large deflections in thin plates and following Simonsen’s 
assumption: 


1 ( dwY 

£. =- 

2\dr) 


( 11 ) 


Consequently, 



The generalized membrane stress then becomes: 


( 12 ) 


N = ot = 


2 " 


The governing equation is then: 


{Nrw')= 


rCA "•' ) 

2 " 


A2n+1 \ 


= 0 


(«) 


(14) 


Subject to the boundary conditions, this equation has the general solution: 

2» 

w = c,r 2 " +1 + c 2 


( 15 ) 


The constants c, and c 2 can be found be applying the boundary condition 
r=rj. The expression for w then becomes: 

2n 2 n \ 

2)1+] _ y 2/7 + 1 

V ) 


2n + \ _L „ 

w =- r. 2 )>+i tan 6 

In 1 


dw 

dr 


-tan 0 at 


(16) 


where £ is the radius of the expanding plastically deformed 2one. 



m 


The kinematic boundary condition (K.B.C.) is then applied at both r = £ and r = r,. 

v r, dw d£ , n , 

K : ^jT + ~dT W =0 at r=r > (note tv = -tan 0= constant) (17) 

where T is time and r, can be considered a time-like parameter. 

Using this kinematic condition, £ can now be expressed in terms of r,: 


t = r, 1 + 


2/i+l 

2n 2« 


2« + l 


(18) 


As a result, (16) becomes 
2n + l 


w = 


2 n 


-r, tan# 


4« + l 

2/1 + 1 


f \ 


\ r \ J 


ln\ 

2/i+) 


(19) 


The vertical punch force, P, is given as the following: 

P = 2?r, JVtan 9 = -?' |C ‘ > ^ tan *£— 

2 " 


The punch displacement is: 
£ = + r, tan # = 2r, tan # 


( 20 ) 




( 21 ) 


Substituting (20) into (21), the force-displacement relation is developed. 


p _ ^C 0 /(tang) 2w J 

2 " 

As an example, consider the following values: 

C 0 =460 MPa 0 =30° n = .22 / =.045” 


P=48.5kN 


( 22 ) 
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Force (kN) 


The analysis provides the force-displacement approximation show below in Figure 8. 


Approximate Load vs. Displacement Curve 

for Varying Thickness Flat Plates using 120 °Cone 
80 -- 



0 10 20 30 1 40 50 61 


Displacement (mm) 

Figure 8: Analytical Approximation for Conical Dishing with Starter Hole 







Radial Cracking and Petalling 

The theory used in this section is derived in detail by Wierzbicki [1999], He proposes that 
the total petalling work is due to crack propagation, petal bending, and membrane 
deformation. Further, these quantities are shown to be interdependent. 



Figure 9: Theoretical Petalling Geometry 


m 


» 


m 


Due to physical limitations of the testing apparatus, fully dished plates could not be 
pet ailed during the senes of experiments. Consequently, a flat plate was used to approximate 
the deformation and fracture conditions seen during petalling. Figure 6 shows an example 
of a six-petal geometry (n pelal =6). From the figure, the central angle, 0, is given as: 



^petals 


(23) 


and the instantaneous crack length, a, is related to the petal length, 1, by 

/ 

a =- 

cos# 


(24) 


m 
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The material is assumed to be rigid-plastic with an average flow stress, o Q . The fully 


plastic bending moment per unit length, M 0 , is calculated as: 



Ultimately, the rate of bending work in a single petal can be computed as 

W b = 4M—/tan# 

P 

where, 

M = rjM 0 


( 25 ) 


(26) 


(27) 


M is the amplified bending Moment, which accounts for the increased bending resistance of 
a curved plate. In the flat plate approximation, rf= 1. The per petal membrane work is then 
shown to be a function of crack tip opening displacement, 8 t : 


^ _ 3.84^/V /3 /(sing)- 473 
tcosQ 

where p is the instantaneous bending radius of the flaps. 
Adding (26) and (28) to get total work: 

w,=w b + w m 

then, 


(28) 


(29) 


dW _g 
dT dT 


4 A/7 tan ( 3>.UM 0 8' n p 2n {smS)' 

p /COS# 


4/3 


(30) 


The total work is then obtained by integrating the above expression with respect to 1: 


W. = 


_ 4A/7 2 tan# | 3.84M < ,7^ 1/ V 2/3 (sin ^) 
p t cos<9 


- 4/3 


(31) 




Load vs. Displacement for .035" Flat Plate 



Figure 10: Approximate Theoretical Load Displacement Curve for Petalled Plate 

Punch displacements, the measured radius of curvature, the experimental crack tip 
opening distance were all used in estimating punch force. Although a closed form expression 
for P vs. 1 was not developed for the load condition, forces were estimated using changes in 
total work over small increments of displacement. 


2 

4M1 7 tan# 

dS 

P 


\- 4/3 


t COS# 


(32) 


(Detailed computations can be reviewed in Appendix G.) In this case a piece-wise force- 
displacement curve for each petalled plate was generated using Wierzbicki’s expressions and 
experimentally measured inputs. An example of a generated curve is shown in Figure 10. 





Strain Field Development 









!• 




Figure 12: Incremental Strain Field Rotation 

As shown in Figure 12 above, the strain field can then be rotated into alignment with the 
direction of crack growth by multiplying by the appropriate direction cosine. Although 


principal stresses are not know, the stress ratios (specifically, which governs necking) 

a 22 

can be determined from e n and e 22 using the Von Mises yield condition and the associated 
flow rule. 


d£ u 

= £n_ 

1 cr 22 

de p 

<7 

2 a 

de n 

_ 1 

a 22 

de p 

2 a 

0 


Divide (36) by (37) and solve for — in terms of to obtain: 

°22 de 22 

dfu + ] 

^ii _ ^£22 2 

a 22 ] + i^£lL 

2 de 22 


(36) 

(37) 


(38) 
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Numerical Modeling 

Preliminary Dishing Model in ABAQUS 

ABAQUS was used to model the preliminary plate’s loading condition. For convenience, 
a pressure loading condition (vice a displacement boundary condition) was selected as a first 
approximation. In this model, a uniform pressure was first distributed over the entire plate. 
This loading geometry resulted in dishing without reverse curvature and a necking fracture at 
the plate s clamped boundary. The uniform pressure, full radius loading analysis was not 
pursued further, since the load condition result in a realistic deformation mode. 


Radius (mm) 



Figure 13: Displacement vs. Radius for Half-Radius Pressure Loading 









Figures 14 and 15 show the displacements for the half radius pressure distribution. It 
was found that using a pressure distribution over half the plate’s radius more closely 
modeled the conditions of hemispherical punch dishing. 






Figure 14: Sectional View of Deformed Plate at 240 MPa Pressure Amplitude 



Figure 15: View of Deformed Plate at 240 MPa Pressure Amplitude 










LS-DYNA Model 


Towards the end of testing, a Finite Element model was created in LS-DYNA using the 
displacement boundary conditions of a rigid sphere, shown in Figure 16. Though the 
revised model moved further from the conditions of gas pressure loading (as in the 
preliminary simulation), it more accurately represented the quasi-static experiment. Three 
conditions were modeled: (1) the 75 mm radius sphere/ .055” plate, (2) 50 mm radius 
sphere/,045” plate, and (3) the 20 mm radius sphere/ .035” plate. Computational results 
were then compared with both the approximate analytical predictions (from Simonsen) and 
the experimental values. Data from this analysis can be found in the “Results and 
Discussion” section of this report. A sample LS-DYNA input file for theR=75mm/.055” 
plate test case can be found in Appendix H. 
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Figure 16: LS-DYNA Simulation with Rigid Body Sphere Boundary Condition 



Design and Testing 


Apparatus Design 

A specialized fixture was fabricated for indentation testing of thin square sheets. 
Designed to work using a 200 KN load cell in a Universal Testing Machine, the fixture 
(shown below in Figure 17) can safely accept a central point load (P) up to the 200 KN 
machine limit. Edge fixity is achieved through the use of serrated bolting ring mating 
surfaces. Detailed size and material specifications for the fixture are shown in Figure 18. 
Further design explanations can be found in Appendix D. 



Figure 17: Notional Test Fixture and Completed Design 

To create the dishing effect, three hemispherical indenters were used with radii of 20 mm, 
50 mm, and 75 mm. Indenter speed (V) was set at a baseline of 10 mm/min. Indentation 
was continued through circumferential cracking until the instance at which radial fracture 
began. At that point, the dishing portion of the test was stopped. A 120° conical indenter 
was then used to simulate petalling from the flat plate condition for the same thickness. 
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Specimen Design 



Figure 19: Test Plate 

Figure 19 shows a representative test plate. Three plate thickness were tested (.035”, 
.045”, and .055”). The plate material has an average uniaxial tensile yield strength (a y ) of 
29,000 psi (200MPa), an average ultimate tensile strength (aj of 46,000 psi (317 MPa), and a 
Rockwell Superficial Hardness of 72. Detailed matenal characteristics can be reviewed in 
Appendix E. The uncoated side of the plate is scribed with a 2.5 mm strain mesh (slightly 
larger than 2 plate thicknesses). During deformation, the horizontal and vertical local 
incremental strains are directly recorded using digital photography. 



36 




















































Testing 



Table 1 shows the test matrix of the series of experiments performed with the Universal 
Testing Machine. Clearly, there were many testing variables to consider. With a finite 
number of test specimens, however, it was a goal of the study to develop a testing 
methodology and identify the most significant factors. 

Tests 1 through 3 were spherical indentation tests, which were used to verify Simonsen’s 
predictions. Tests 4 through 7 were used to determine the best approach given the 
limitations on time and equipment. In Test 4, a conical punch was used to indent a plate 
with a preformed hole. Later, Tests 14 andl5 were used to observe to effect of hole size in 
this same loading condition. Test 5 was the first experiment in which the starter hole was 
pre-cracked in order to promote the petalling behavior. Tests 8, 9, and 10 also used pre¬ 
cracked holes to examine the effect of plate thickness. Tests 11 and 12 investigated the 
effect of increasing the number of pre-cracks, while Test 13 was used to observe the effect 
of increasing pre-crack length. Tests 4,14 and 15 confirmed the effect of starter hole size. 
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Results and Discussion 


Spherical Indentation 

Figure 19 shows the experimental and approximate analytical results for spherical 
indentation testing. (Note, the analytical approximation is valid only up to Simon sen's failure 
prediction. Points beyond the maximum load are artifacts of the computation.) 


Spherical Indentation Results 



Cross Head Displacement (mm) 


Figure 20: Approximate and Experimental Results for Spherical Indentation 

Figures 21 through 23 on the following page, show results of the LS-DYNA numerical 
simulations for the same three load cases. 








Spherical Indentation Results for R=20mm Sphere / .035" Plate 
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Figure 21: LS-DYNA Results for R=20 mm punch /.035” Plate 


Spherical Indentation Results for R=50 mm Sphere / .045" Plate 
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Figure 22: LS-DYNA Results for R=50 mm punch /.045” Plate 


Spherical (dentation Results for R-75 mm Sphere / .055" Plate 
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Figure 23: LS-DYNA Results for R=75 mm punch /.055” Plate 
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As in Simonsen’s experiment, the experiments correlated to the analytical approximation 
within 5%. The preliminary LS-DYNA model (pictured in Figure 24) also resulted in good 
agreement, although the model requires further refinement for follow-on tests. This 
portion of testing validated the testing method and clearly proved that the new specimen test 
fixture functioned as designed. 

During spherical indentation testing, it became clear that the specimen material was too 
ductile to carry out all tests as planned. Plastic deformations were quite large during 
dishing, and the shallow conical punches did not provide sufficient depth of stroke to 
continue petalling. Furthermore, the dishing failure always resulted in incomplete 
circumferential cracking (as in Figure 24 below). As a result of these two obstacles, the 
dishing and petalling phases were modeled in two entirely separate experiments. 


■m 


150 mm Diameter/.055” Plate 


40 mm Diameter/.035” Plate 






Figure 24: Experimental and LS-DYNA Spherical Dishing Deformations 
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Conical Indentation without Pre-cracks 


Figure 25 shows the force-displacement relation observed during conical dishing. During 
this experiment, a constant plate thickness was tested while the starter hole was varied in 
size. The dashed line represents the analytical approximation derived earlier. The analytical 
approximation suggests that the force-displacement curve should be independent of hole 
size. This behavior was observed up to a displacement of about 20 mm. After that point, 
the 3 experiments diverge. Certainly the smaller hole size induced much earlier radial 
cracking due to much higher hoop stresses. 

It is difficult to draw immediate conclusions, since each case’s geometry, and 
consequently the material’s load history, cannot be directly correlated. The analytical 
approximation, however, serves as a broad estimate of the material’s general force- 
displacement behavior, but falls short in accounting for the problem’s changing geometry. 


.035" Plate Dished with 120° Cone 

(10, 22, and 45 mm Hole, No Precracks) 



Figure 25: Conical Indentation Force-Displacement Curves (with Starter Hole) 








Conical Indentation With Pre-cracks (Petalling) 


Load vs. Displacement for Varying Thickness Flat Plate 
120° Cone, Four 3-mm Pre-cracks, 9.5 mm Hole 


Initial Condition: 9.5 mm Hole with 4 Pre-cracks 



Displacement (mm) 
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Figure 26: Conical Indentation Force-Displacement Curves (with Pre-cracks) 


Figure 26 shows the results of the plate petalling experiment, which was conducted with a 
9.5 mm starter hole and four, 3 mm pre-cracks. Figure 27 shows how the petalling test 
progressed. 



10 mm 



25 mm 



55 mm (completed test) 


Figure 27: Petalling Test Punch Displacement Progression 






The petalling test showed that the effect of plate thickness was approximately linear. In 
each experiment during the first 10 mm of displacement, the cone settled into the starter 
hole and propagated the pre-cracks. After that point, the petals began to form. From the 
pictures, one can see that very little curvature developed. Again, this was due to the shallow 
angle of the conical punches. Only 55 mm of punch displacement could be observed before 
the punch began to interfere with the fixture’s boundary conditions, at which point the test 
had to be stopped. 

Although, there appears to be disagreement between the analytical approximation and the 
observed force-displacement curves, the general trend and level of force correlate reasonably 
well. It is apparent that the artificiality of pre-cracking has some impart on the experiment. 
In the two experiments where natural radial cracks formed and petalling could be observed 
(Figure 28), the force-displacement curve closely resembled the behavior derived out of 
Wierzbicki’s expressions Eq.(32). 



Displacement (mm) 



Pre-Dished and Cut Specimen Conical Dishing with Starter Hole 

Figure 28: Force-Displacement Observations During Petalling 








It is important to keep in mind that the dishing force is five to six times higher than the 
less critical petalling force. With the perspective that this approach is simply another 
available tool, this petalling approximation could certainly provide useful first-order 
estimates to a designer or post-explosion investigator. 

Effect of Pre-crack Length 

Figure 29 shows the effect of pre-crack length. Variation in the length of the starter 
crack produced little effect, other than lengthening the amount of time required for the cone 
to settle into full contact with the plate. The result was a small shift in the force - 
displacement curve. As displacement increased the effect was diminished. This behavior 
was expected, but required verification. The length of pre-crack, however, did affect the 
number of petals formed. The 3-mm pre-crack yielded four petals, while the 10-mm pre¬ 
crack yielded only four. This behavior agrees with Wierzbicki’s prediction than some energy 
absorption minimum exists, tending to promote the formation of three to five petals. 
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Effect of Pre-Crack Length 
Eight Pre-cracks, .035" Plate, 9.5 mm Hole 



Figure 29: Effect of Pre-Crack Length on Force-Displacement Curve 













Effect of Number of Pre-cracks 

As shown in Figure 30, the number of pre-cracks had no significant effect of the overall 
force-displacement curve. Only four petals developed in each case, however. In both tests, 
all pre-cracks tended to start, but only four continued to grow. As in the previous 
examination of pre-crack length, this result also suggests that a petalling energy minimum 
exists, which is independent of the number of initial radial necks or cracks. 


Effect of Number of Precracks 
3 mm Pre-cracks, .035" Plate, 9.5 mm Starter Hole 



Cross Head Displacement (mm) 


Figure 30: Effect of Number of Pre-Cracks on Force-Displacement Curve 
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Crack Displacement and Strain Field Examination 


Figure 31 shows the increment of crack growth used while mapping the strain field 
around the crack tip. The crack was allowed to grow approximately 2.5 mm between framed. 
The white arrow gives a fixed point from which to visually reference the crack’s advance. 

The mapping grid was laid out according to the method discussed in an earlier section. 



Frame 1 Frame 2 

Figure 31: Crack Advance Increment 


The displacement filed vectors are plotted relative to the crack tip. Examination of the 
displacement field shows that, at this instant in time, the membrane is rapidly stretching and 
thinning in the region in front of the crack. 

Biaxial incremental strain values were computed for a small rectangular patch around the 

crack tip. Similarly, the principal stress ratio, —, was computed from Eq.(38) for each of. 

<J 22 

the grid points, helping to identify regions where necking becomes critical. Figure 32 shows 
the measured displacement field in the postage stamp size area surrounding the crack tip. 
Figure 33 shows the measured incremental strain fields and Figure 34 shows the computed 




field. A visual inspection of the displacement and strain fields suggests that necking 

fj 22 

occurs primarily directly ahead of the crack. This behavior was both expected and observed 
during the experiment. 
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Figure 32: Crack Advance Incremental Displacement Field 
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Figure 33: Resultant Strain Increment Fields 
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Conclusions and Recommendations 


Conclusions 

This study has investigated using simple, inexpensive quasi-static tearing tests to gain 
insight into the complex dynamic problem of explosions. It has been shown that these tests 
can successfully reproduce plate deformations similar in nature to those found after large- 
scale explosions. It was also shown that the forces, displacements, and energy dissipation in 
the plate-dishing phase can be accurately approximated by Simonsen’s approach (within 5%). 
Although this series of experiments was not able to generate spiraled petals, Wierzbicki’s 
calculation of petalling work appeared to provide reasonably good (within 20%) estimates of 
energy dissipation during the second phase of deformation. The conical punch test, 
however, requires a revised mechanical design in order to better model the flat plate petalling 
and post-dish petalling cases. 

An approach to developing fracture criterion was also presented by measuring 
incremental strains around the crack tip. Subsequently, the local stress ratios governing 
material necking were computed. With that information in hand, crack growth can be 
predicted in the next increment. This approach could result in a revised crack growth 
criterion, which might then be used to improve existing hydrodynamically coupled finite 
element simulations. 
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Recommendations 


• Th e punch tests should explore a wider range of penetrator geometries. Punch variation 
may introduce more curvature into the resulting petals, helping to capture more of the 
dynamic effects. Pyramid or volute shapes might be considered. Certainly, narrower 
cone angles should be tried in order to increase the punch displacement and extend 
observation time. 

• Future analytical work should attempt to more closely link the quasi-static test results 
with characteristics of the dynamic problem. To do so may require dynamic testing, 
which could be accomplished at government blast testing facilities. 

I 

• Future numerical simulations should apply new punch geometries and include new plate 
structures (sandwiches, foam, etc...) which might potentially be considered for ship hull 
reinforcement. Future numerical simulations should also attempt to more completely 
capture the dynamic aspects of an explosion. 

• The potential improvement in crack growth criterion could result in a tremendous 
advance m the understanding of crack propagation due to intense blast pressures. This 
avenue should be explored more fully in order to update current prediction codes. 

• Other areas of application for this research should be investigated. Certainty pipeline 
explosions (as in the natural gas industry), internal aircraft explosions, and pressure 
vessel explosions are just of few items for consideration 
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Appendices 




Appendix A: Introductory Material 


Since the inception of the warship, designers have sought means to make combatant 
vessels more damage resistant. Only since World War II, however, has there been significant 
effort to compile and analyze detailed reports of vessel battle damage. In the subsequent 
post war years, much analysis was conducted in the areas of hull armor, resistance to mine 
attack, and torpedo side-protection. As might be expected, the systems that were ultimately 
developed to enhance survivability, incurred substantial weight, space, and performance 
penalties. As the years have passed, the damage database has deteriorated. Since the late 
1940’s, the U.S Navy has relied largely on full-scale tests or “SINKEX’s” to demonstrate 
battle damage effects. There have been only a handful of wartime incidents demonstrating 
the realistic damage absorbing and recovery capabilities of modem warships. In that time, 
threat weapons (mines, missiles, and torpedoes) have been improved markedly. In contrast 
to their World War II predecessors, state of the art weapons are extremely likely carry out a 
successful attack and inflict extraordinary damage with a high degree of accuracy. 

As seen more recently in the case of the U.S.S. Cole, tremendous hull damage can be 
suffered during unconventional close aboard explosions. Figure A1 shows the 20-ft by 40-ft 
amidships hole torn in the ship’s port side by a waterline explosion. 

Inspection of the picture shows that the explosion created a spherical bulge, or dishing, in 
the ship’ hull, prior to tearing into numerous petals. Though officially unconfirmed, it has 
been suggested that the attack was carried out using between 400 and 700 lbs. of C-4 military 
explosive. Unofficial standoff distance estimates range from on-contact to 10 feet. 




Figure Al: USS Cole Port Side Damage 


In this instance, a majority of the damage occurred below the waterline. Figure 2 clearly 
shows the above-waterline tearing. From Figure A2, the characteristic petal formations car 
be identified, as well as the upper deck hard point at which crack propagation is arrested. 
Hull panels at this location on a typical naval combatant ship are usually constructed of 
either mild steel (MTS-45) or HY-80, and they can vary between .375” and .5” in thickness 
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Appendix B: Additional Background 

Holing prediction for metal plate has been a topic of study since the early 1900’s. The 
research is rooted in a 1912 study done by Bertram Hopkinson, which examined the 
resistance of naval plating against artilleiy shell penetration. In that case, the deformed 



Figure Bl: Hopkinson's 1912 Armor Plate 

plating, shown in Figure Bl, developed a coronet-shaped pattern with radial cracks and 
necking at the petal edges. Since that time there have been countless experiments 
investigating armor effectiveness. Until the past two decades, the vast majority of plate 
holing and deformation studies have focused on the impact and associated damage of 
ballistic particles. The most notable and comprehensive exception is the 1940’s work of Sir 
G.I Taylor (1948) who examined the nature of submerged blast waves and their effects on 
thin plates. Taylor’s work in the area of submerged explosions was continued by Robert 
Cole (1948) in the text, Underwater Explosions . Together, these two researchers provide 
much of the technical underpinning for current blast damage prediction methods. 
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Explosive deformation and holing studies in naval vessels, largely classified, have been 
empirical derivations based on years of accrued data. An undisclosed solution developed by 
ONR, predicts a minimum and maximum expected hole radius in naval panels for a given 
set of panel parameters. This engineering tool was the result of hundreds of live fire tests, 
and is one of the most commonly used methods used by U.S. Navy designers in estimating 
expected battle damage due to close proximity explosions. This solution suggests the 
following general relationship: 


Stiffened Panel 



Figure B2: Current Panel Damage Formulation 


^Min — — ^Max = /(G, T, M) 

R), = Resultant Hole Size 

R™, = Minimum Predicted Hole Size 

R^ = Maximum Predicted Hole Size 

G = Panel Geometry 

T = Plate Thickness 

M = Material Properties 
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Other problems, however, lend themselves to the study of panels under explosive loading. 
Wierzbicki and Thomas (1993), by examining plate cutting behavior, were able to related the 
geometiy of the deformed plating geometry to the mechanical and frictional work done in a 
vessel grounding, thereby predicting damage. Later, Wierzbicki (1996) and Nurick (1996), 
investigated the response of clamped thin sheets when subjected to on-contact charges. 
Wierzbicki (1999) subsequently proposed that the kinematics of the cutting process, shown 
in Figure B3, were very similar to both those of the explosive petalling and ballistic 
penetrator problem. 



Figure B3: Cutting and Petalling Geometric Similarities 

Along similar lines, Atkins (1998) observed the necking and radial crack formation in 
ductile materials when perforated by both spherical and conical penetrators. Using a 
hydraulic bulger, he later analyzed thin sheet necking around biaxially loaded holes (1999). 
Simonsen (2000), built upon this work, using spherical indenters to analyze force- 
displacement relationship in clamped mild steel plates. Again, these studies were targeted at 
making damage prediction for vessel collisions and ballistic penetrators. Nevertheless, these 





experiments closely model thin sheet behavior during the dishing, disking, and pre-petalling 
phases of an explosion. 

To date, there has been little effort to correlate the results of quasi-static indenter 
experiments to the damage observed both in live fire explosive tests and in combat. With 
that goal in mind, this research is deemed both relevant and unique. 




Appendix C: Problem Statement Details 


Currently, U.S. Navy surface ships are built using a damaged length design standard 
equivalent to 15% Length Between Perpendiculars (LBP). This rule for damaged length is 
based on both past operating experience and on the damage expected from a given threat 
weapon. As illustrated in Figure Cl, accurately determining the opening size is crucial in 
determining compartment layout, equipment placement, and bulkhead spacing. In a poorly 
designed 500-ft ship, for instance, a 15% damage length (75ft) could result in total loss of 
propulsion. By the same token, over designing a ship using excessive damage length results 
in increased structural weight, higher cost, and decreased performance. 



Figure Cl: Surface Ship 15% LBP Damage Length 

Damage prediction methods can be extraordinarily complex. The most modem 
techniques use machine codes, which couple explosive mechanics with the concurrent 
hydrodynamic effects. These codes, however, are both expensive and time consuming. 
Furthermore, the accuracy of their output is not yet fully validated. Other than scaled live 







fire tests, no reliable prediction tool exists which can quickly and accurately estimate damage. 
An analytical method which uses inputs such as charge size, standoff distance, material could 
prove useful both to ship and weapon designers. 




Appendix D: Fixture Design Supplemental 


A significant amount of design work was required to ensure that the mechanism (1) did 
not fail under the expected test loads and (2) was adaptable for use in the Universal Testing 
Machine (UTM) belonging to the MIT Impact and Crashworthiness Laboratory. 

As shown in the drawings, the fixture horizontally holds a 300 mm x 300 mm square thin 
sheet while an indenter or punch is lowered. The clamped dimensions, shown in Figure Dl, 
are 220 mm x 220 mm. The specimens maximum dimensions were selected based on 
several criteria. First, the specimen needed to be large enough to easily observe the resultant 
deformations, yet small enough to fit within the UTM. Further, it was desirable to have a 
portable, light-weight test fixture (in this case 75 lbs. was designated as an upper weight 
limit). Lastly, steel sheet stock is normally manufactured in 4 ft x 8 ft sheets. A 300 mm 
square sample allows 32 specimens to be taken from a single sheet with minimal waste. 


220 mm Clamped Width 


I 


1.5 mm Fillet Radius 


300 mm Specimen Width 
(1.15 mm Specimen Thickness) 


10 mm 


25 mm 


Figure Dl: Clamped Specimen Cross Section 




Since the test specimen represents a single clamped panel, the specimen needed to be 
geometrically similar to existing naval panels. Frame spacing within a mid-sized combatant 
is generally near 8 ft, with the panels being nearly square (or slightly rectangular). At 
thicknesses of 1/2” to 3/8”, a panel width-to-thickness ratio (w/t) of 200 can be reasonably 
assumed. In order to maintain this ratio, three specimen thicknesses (.89mm, 1.14 mm, and 
1.40 mm ) were selected from available mild steel stock. 

The body of the fixture is machined from mild steel (45ksi Yield Strength). Member 
thicknesses were selected to prevent yield during normal testing and to minimize weight. 
Note that the fixture narrows to 6” at its base to transfer force directly to the load bearing 
beam in the UIM. Thin lateral bracing panels were used to prevent fixture racking during 
testing. The fixture mounts to the UTM through a series of threaded holes in its baseplate. 

The upper and lower bolting rings are used to hold the specimen in place. These rings are 
machined from A2 Air Hardened Steel (108 ksi Yield Strength). This grade of steel was 
selected to satisfy fixity requirements. Firm edge fixity is required in order to make 
boundaiy condition assumptions. Past experiment have shown that simple bolting flanges 
do not provide sufficient clamping restraint. Typically a drawbead can be used to overcome 
this problem. When testing thick or high strength materials, however, it becomes difficult to 

j 

form the specimen around a drawbead. For this reason (and since the fixture is intended to 
test a wide range of materials and thickness), a serrated mating surface was used to provide 
fixity. 

In order to ensure sufficient tooth hardness and clamping force, an estimate of material 
strength was made using Slip Line Fracture Mechanics. The analysis showed that the mating 
surfaces must have a yield strength at least 50% higher than that of the test specimen. The 
estimates are summarized in Figures D2 and D3. 




Test Material 

. X 



,45 Degrees 


mm 


Figure D2: Mohr Circle For Tooth Element Tooth 


■cr rf =-tf-k=-2k ( (l+^) 


tf r W = -2k w (l + f) 


d = -k-2k 


Figure D3: Mohr Circle For 45-Degree Rotation Tooth Element 


= ^n, 


M ,+ f> = Mi+f) 

k ] + ~ 

T~ = “■ = 1 -44 (Ratio of Base Plate Yield to Test Material Yield) 





Estimate of Frame Sizes and Stresses Induced During 
Indentation Testing 


• Cross Beam Dimensions: 

h = 25.4mm b := 50-mm L:= 300mm 

Expected Load and Material Characteristics: 

P := 10500CN P = 2.36 x 10 4 lbf 

Y := 108000psi Y = 7.446x 10 8 Pa 


Assuming a Rigid Test Specimen and a 

Simply Supported Beam with a Concentated Central Load: 


M := 


L P 
2 4 



a beam ‘ j a beam _ ^.324x 10 Pa 

Safety Factor: 

Y 

-= 1.017 

a beam 

...or distibuting load along ther beam, gives the same result for the maximum 
developed stress: 

i 

P ^ 1 T 2 

w:=— M := — • vv-L 

4L 2 

°beam‘ = ^’Y c beam~^'^^ x ^ Pa 

Safety Factor: 

Y 

-= 1.017 

a beam 
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Using the exact solution for sinusoidally loaded, simply supported plate: 

F D ; =-7 M:=F D i 

71 

CT beam := M • — 

Safety Factor: 

Y 

-= 2.509 

a beam 

Compressive Load on Legs: 

Yjeg := 45000 psi 

w , eg := 49.3 mm A leg : = 2 ( 2 t|eg- w Ieg) 

<T| eg = 2.662 x 10 7 Pa 



CT beam — 2-968 x 10 Pa 
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Appendix E: Uniaxial Testing 


The mild steel plate material was tested Uni.-axially in tension according to ASTM 
standards. The .035” and .045” were tested in 1999 in concert with a separate project. The 
.055” thickness was tested specifically for this research. All specimens were tested in the X 
and Y direction to verify that the material behaved isotropically. 


Load vs. Time for .035" X Sample 2 



125 175 225 275 325 375 425 


Time (Sec) 

Figure El: Representative Load vs. Time Curve 


Load Change vs Time .035 X Sample 2 



Figure E2: Representative Load Change vs. Time Curve 





During the uni-axial testing, all of the thin strip specimens displayed stable fracture, 
consequently the precise load at fracture was not readily apparent. By plotting the change in 
load vs. time for each specimen, a “knuckle” in the load vs. Time plot was clearly identified, 
which indicated the onset of fracture. Figures El and E2 illustrate this approach. Figures 
E3 and E4, respectively, show a representative uni-axial specimen and a close-up view of the 
fracture. 


















The final dimensions of each specimen’s cross section were measured in order to 
determine the true stress and strain of the specimen at fracture. This data point was vital in 
constructing the true stress-strain curve for the finite element simulation. Due to necking,. 
the specimen’s highest true fracture stress and strain were achieved in the middle, where the 
fracture initiated. Figure E5 shows the approximate final geometry for the necked 
specimens. 



Figure E5: Approximate Necked Specimen Geometry 


C* 

Reduction in Squared Thickness (RST) and true fracture strain for the thin strips were 
computed in accordance with McClintock and Zheng (1993), 


RST = \-{t f lt 0 f 
s s = -ln(\-RST) 


(El) 

(E2) 
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m 


dtf 



Figure E6: Representative True Stress vs. Strain Curve for Uniaxial Test 


Some formulation was required to compute o nilM[ (o^J : 


£e)~r 

as 


(E3) 




Pfracture (E4) 

Where substituting (E3) into (E4) yields, 

^'fracture ~ ^m^ni 4" (®"m — (^m ~ ^e)~l ) 


m 

(E5) 


and solving (E5) for o m 

rr ___ as 


A, + A 


(E6) 






Table El: Edge and Middle Thickness at Fracture for Uniaxial Tests 
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Table E2: Right and Left Edge Fracture Strain for Uniaxial Tests 


Sample 


.035 X#1 


.035 X #2 


Right Edge 

V Right Edge Fracture Strain 






Fracture Area (mm^) 







.055 Y #1 


.055 Y #2 

0.7975 | 


Sample 


.035 X #1 


.035 X #2 



Left Edge 

Left Edge Fracture Strain 


1.5995 


1.5019 




.035 Y #1 


.035 Y #2 


0.7980 


0.7666 


1.5995 


1.4549 


Fracture Area (mm*) 


1.1493 


1.2180 


1.1200 


1.1768 


























































































































Table E3: Middle Fracture Strain for Uniaxial Tests 



# 
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Table E4: Estimation of Average Uniaxial Fracture Stress 




Computed Values 


■ ; Av 


Sample 


Edge Fracture Strain (Epsilon E ) 


.035 X#1 


.035 X #2 


1.705 


1.602 


Total Edge Area (A E ) 


2.1837 


2.3200 


dSigma/dEpsilon 


300 


300 


.035 Y #1 








1.625 


2.2120 


300 


.035 Y #2 


1.478 




2.3262 


300 






.045 X #1 


1.422 


2.9307 


300 


.045 X #2 


1.369 


3.0782 


300 








.045 Y#1 


1.440 


2.6788 


300 


.045 Y #2 


1.458 


2.7353 


300 






.055 X #1 


1.582 


2.8829 


300 


.055 X #2 


1.550 


2.9475 


300 


.055 Y #1 








1.597 


2.8896 


300 


.055 Y #2 


1.613 


2.8646 


300 


mmm 








Sample 


Fracture Load (lb) 


Average Fracture stress (MPA) 


Fracture Stress Middle (MPA) 


.035 X #1 


697.78 


1019.08 


1120.18 


.035 X #2 


696.04 


953.24 


1048.34 








































































































Table E5: .035” Uniaxial Sample Test Data 


.035" Thickness X Direction 


file name: c2x-1.dat 
description: 0 035" sheet metal, "X" direction 
'"' test date: 11/9/1999 
test mach.; tnstron 1331 
specimen #: C2X-1 
area: 0.0176 
thickness: 0.0350 
.width: 0.5035 


file name: c2x-2.dat 

description: 0 035" sheet metal. "X" direction 
test date 11/9/1999 
test mach.; tnstron 1331 
specimen #: C2X-2 
area: 0.0176 
thickness: 0 0350 
width: 0.5030 


pretest notes 

data acquisition - 300 msec, 22 bits 
channel 1 = strain: 10%/v 
channel 2= load: 562 Ibs/v 

sampling rate = 0.5 sec ._ 

load rate = 0.2 in/min 

first data point = zero load, second data point = zero strain 


C2X-1 

test notes 

yeild stress: 

23,300 

ult stress: 

46,294 

break point: 

inside qaqe area 


.035 Thickness 


pretest notes 

data acquisition = 300 msec, 22 bits 1 
channel 1 = strain: 10%/v 
channel 2 = load: 562 IbsAr 
sampling rate = 0.5 sec 

load rate = 0 2 in/min ! I 

first data point = zero load, second data point = zero strain 


C2X-2 

test notes 

yeild stress: 

23,300 

ult. stress: 

45,903 

break point: 

inside gage area 


Y Direction 


filename: c2y-1 dat 

description: 0 035" sheet metal. "Y" direction 
test date: 11/9/1999 

test mach.; Instron 1331 

specimen#: C2Y-1 

area: 0 0176 

thickness: 0 0350 
width: 0.5040 


file name: c2y-2.dat 

description: 0.035" sheet metal. "Y" direction 
test date: 11/9/1999 

test mach.; Instron 1331 
specimen #:' C2Y-2 

area: 0.0177 

thickness: 0 0350 
width: 0.5045 


pretest notes j. 

data acquisition = 300 
jchannel 1 = strain: 10%/ 
channel 2 - load: 562 


sampling rate =0.5 sec 
load rate = 0.2 in/min 


msec, 

i/v 


. 22 bits 


first data point = 

zero load, second 

data point = zero strain 

C2Y-1 

test notes 

1 

yeild stress: 

25,100 


ult. stress: 

45,176 


break point 

inside qaqe area 



retest notes 


data acquisition = 300 msec, 
channel 1 = strain: 10%Ar 
channel 2 = load: 962 ibs/v 


sampling rate = 0.5 sec 
load rate = 0.2 in/min 

first data point = zero ioad, second data point = zero strain 


C2Y^ 

test notes 

yeild stress: 

25.100 

ult stress: 

45,146 

break point: 

inside qaqe area 
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Engineering Stress (psi) Engineering Stress (psi) 


0.035 Sheet Metal, M X" Direction 



0 5 10 15 20 25 30 35 40 45 * 

% Strain 

Figure E7: Engineering Stress vs. Strain Curve for .035” X Uniaxial Test 


0.035” Sheet Metal, "Y" Direction 



5 10 15 20 25 30 35 40 - 45 

% Strain 


Figure E8: Engineering Stress vs. Strain Curve for .035” Y Uniaxial Test 








Table E6: .045” Uniaxial Sample Test Data 


.045“ Thickness X Direction 


file name: c3x-1 dat 

description: 0 045' sheet metal, "X’ direction 
test date: 1 iy9/1999 

test mach.; Instron 1331 
specimen P: C3X-1 

area: 0.0227 

thickness: 0 0450 
width: 0.5050 


file name: c3x-2.dat 

description: 0 045' sheet metal, “X" direction 
test date: 1179/1999 

test mach.; Instron 1331 
specimen #: C3X-2 

area: 0.0227 

thickness: 0.0450 
width: 05050 


pretest notes 

data acquisition = 300 msec. 22 bits 
channel 1 = strain: 10%/v 
channel 2 = load: 562 Ibs/v 
sampling rate = 0 5 sec 
load rate = 0.2 in/min 

first data point = zero load, second data point = zero strain 


pretest notes 

data acquisition = 300 msec, 22 bits . 
channel 1 = strain: 10%/v 
channel 2= load: 562 Ibs/v 
sampling rate = 0.5 sec 
load rate = 0.2 in/min 

first data point = zero load, second data point = zero sir; 


OX-1 

test notes 

yeild stress: 

28,600 

ult. stress: 

47,776 

break point: 

inside gage area 


C3X2 

test notes 

yeild stress: 

28,600 

utt stress: 

47.688 

break point: 

inside gage area 


.045" Thickness V Direction 


file name: c3y-1.dat 

description: 0.045" sheet metal, 'Y“ direction 
test date: 11/9/1999 
test mach.; Instron 1331 
specimen #: C3Y-1 
area: 0.0226 
thickness: 0.0450 
width: 0.5030 


file name: c3y-2.dat 

description: 0 045" sheet metal. "Y" direction 
test date: 11/9/1999 
test mach.; Instron 1331 
specimen #: C3Y-2 
area: 0.0226 
thickness: 0.0450 
width: 0.5030 


pretest notes 

data acquisition = 300 msec, 22 bits 
channel 1 = strain: 10%/v 
channel 2 = load: 562 Ibs/v 
sampling rate = 0 5 sec 
load rate = 0.2 in/min 

first data point = zero load, second data point = zero strain 


C3Y-1 

test notes 

yeild stress: 

30,000 

ult. stress: 

46,948 

break point: 

inside qaqe area 


pretest notes 

data acquisition - 300 msec. 22 bits; 
channel 1 = strain: 10%/v 
channel 2 = load: 562 Ibs/v 
sampling rate - 0.5 sec 
load rate = 0 2 in/min 

first data point = zero load, second data point = zero 


C3Y-2 

test notes 

yeild stress: 

30,000 

ult. stress: 

46,999 

break point: 

inside gage area 
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Engineering Stress (psi) 


0.045” Sheet Metal, "X" Direction 



Figure E9: Engineering Stress vs. Strain Curve for .045” X Uniaxial Test 


0.045" Sheet Metal, "Y" Direction 








Table E7: .055” Uniaxial Sample Test Data 


file name: c4x-1 dal 

description: 0 055“ sheet metal. X direction 
test date: 6/22/2002 

test mach.; UTS 

specimen #: C3X-1 

area: 0 0278 

thickness: 0 0550 
width: 0.5050 


.055” Thickness X Direction 

file name: c4x-2 dat 

description: 0 055“ sheet metal. X direction 
test date: 6/22/2002 
test mach.; UTS 
specimen #: C3X-2 
area: 0 0278 
thickness: 0 0550 
width: 0 5050 


pretest notes 

data acquisition = 300 msec. 22 bits 
channel 1 = strain: 10%/V 
channel 2 = load: 562 IbsAr 
sampling rate = 0.5 sec 
load rate = 0.2 in/min 

first data point - zero load, second data point - zero strain 


pretest notes 

data acquisition = 300 msec, 22 bits 
channel 1 = si rain: 10%/v 
channel 2 = load. 562 IbsAr 
sampling rate = 0 5 sec 
load rate = 0.2 in/min 

first data point = zero load, second data point = zero strain 


C3X-1 

test notes 

yeild stress: 

8 

oo 

R 

ult. stress: 

48,202 

break point: 

inside gage area 


OX^ 

test notes 

yeild stress: 

28.900 

ult. stress: 

48,321 

break point: 

inside gage area 


.055” Thickness Y Direction 


file name: C4Y-1 dat 

description: 0 055“ sheet metal, *Y“ direction 
test date: 6/22/2002 
test mach.; UTS 
specimen #: C4Y-1 
area: 0.0277 
thickness: 0 0550 
width: 0 5040 


pretest notes 

data acquisition = 300 msec, 22 bits 
channel 1 = strain: 10%Ar 
channel 2 = load: 562 IbsAr 
sampling rate = 0.5 sec 
load rate = 0.2 in/min 

first data point = zero load, second data point = zero strain 


C4Y-1 

test notes 

yeild stress: 

30,tOO 

ult stress: 

47,129 

break point: 

inside gage area 


file name: C4Y-2.dat 

description: 0.055’ sheet metal, "Y“ direction 
test date: 6/22/2002 
test mach.; UTS 
specimen U: C4Y-2 
area: 0 0277 
thickness: 0.0550 
width: 0 5040 

pretest notes , ] 

data acquisition = 300 msec. 22 bits 

channel 1 - strain: 10%/v 

channel 2 = load: 562 lbs/v 

sampling rate = 0.5 sec 

load rate = 0.2 in/min 

first data point = zero load, second data point = zero strain 


C4Y.2 

test notes 

yeild stress: 

30.215 

ult. stress: 

47,111 

break point: 

inside gage area 






























Appendix F: Analytic Solution Worksheet 
.035" Thickness using 20 mm Radius Punch 


Entering Arguments for Analytic Solution 

Work Hardening Exponent and Coefficient: 

n := .22 C o :=46010 6 Pa 


Material Yield and Ultimate Strength 


a y := 23300psi t? u := 46100psi 


Plate Thickness 


t Q := .035in 


Plate Clamped Radius R j \$ mm 


Punch Radius 


= 20 mm 


Compute Failure Parameters 

Vc_fail := <- 957 + -399 n)-rad V c _f ail = 1 045 

4R <7 

%ail := l*41’ n 5^ a jj = 38.781mm 


Compute Stresses 


Flow Stress 



+ °u) 


cr 0 = 2.392x 10 8 Pa 



Generalized Membrane Stress 




Plate Predictions at Failure 


R := 110 


P fai! : = 2 N 0 R b( sin (v c _fail)) 2 _ 


Pfail = 2-308x IO N 


R cfail := R b sin (^c_fail) 


Ref := —=— 

mm 


r<-:= Ref.. R 


w fail( r f) : = - 


Pfaif'^ 


w 


R + 


,j( r f ) 2 - ( sin (^ c fail)) 4 ,] 


j(R ) 2 ~ ( sin ( v t / c_fail )) 4 


2-Tt-N 


mm 


Establish Load Range in 5000 N Increments, Up to the Predicted Failure Load: 
Load: 

i:= 1..7 


P := 


0 > 

p o 

:= ON 




P fail 





5-N 


Pfail 




P 1 

5 

p i = 

= 4.616x 

10 N 

P fail 






4-N 

P 2 

Pfail 

4 

p 2 = 

= 5.771 x 

io 3 n 

P fail 



i 



3-N 

•N P 3 

Pfail 

3 

P 3 = 

= 7.694 x 

io 3 n 

P fail 






2*N 

P 4 

P fail 

P 4 = 

= 1.154x 

io 4 n 

2P fail 


Z 




3-N 

Pfail 

P 5 

_ 2P fail 

3 

P 5 = 

= 1.539x 

io 4 n 


V N / 




Wrapping Angle 

( 0 
22.753 
25.62 

ij; j = | 29.954 deg 
37.699 
44.92 
59.861 ) 


\|/j := asin 


( I p ^ 


2-"-N 0 -Rb 


\|/0 := asin 


vj/2 := asin 


\|/4 := asin 


2' 7t 'N 0 -R b 


2-n-No-Rb 


2-«-N 0 -Rb 


\(/6 := asin 


fail 
1 2-7t-N 0 Rb 


v|/l := asin 


\[/3 := asin 


2*N 0 Rb 


2 ' 7t 'N C) R b 


vg5 := asin 


2’tN 0 .R b 





Strain Components 


E rr := ,n 


e tt '= - £ rr 


z n = 0.143 


-0.103 

-0.143 

-0.234 


-0.689 


Deformed Plate Thickness 


t := t 0 -cosf vpjj 


True Material Stress: 


°rr :=C o(- ,n ( C0S ( v t y i))) n 


2.646x 10 


2.793 x 10 


o n =\ 3x 10 Pa 


3.342 x 10 


3.64 x 10 


4.238x 10 
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Strain Components 


z n := In| 


1 


5 M 


V C0S V ' V J 


rr 


£ tt £ rr 


f 0 > 


f 0 

0.081 


-0.081 

0.103 


-0.103 

0.143 

s tt “ 

-0.143 

0.234 


-0.234 

0.345 


-0.345 



^-0.689; 


Deformed Plate Thickness 


True Material Stress: 


t :=t 0 -cos(\j/j) 

a rr “ C o ‘ ( ln ( cos ( V i))) n 


a rr ~ 


2.646x 10 


,8 


2.793 x 10 


8 


3 x 10 


8 


3.342x 10 


3.64 x 10 


8 


4.238x 10 


.8 


Pa 
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Solve for Shape of plate, valid for 1^. < r < R 

R:= 110 



0 

4.616 

5.771 

7.694 

11.541 

15.388 

23.082 


N 10 3 



0 

0.397 

0.447 

0.523 

0.658 

0.784 

1.045 



0 

22.753 

25.62 

29.954 

37.699 

44.92 

59.861 


deg 



RcO:»R b . Si “ (v0) 

mm 

RcO = 0 


Rcl:=R b " n( ' l ’ l) 

mm 

Rcl = 7.735 

# 

R c 2»R t , sin(v2) 

mm 

Rc2 = 8.648 


Rc 3: .R k . iin< ' f3) 

mm 

Rc3 = 9.986 

# 

Rc4:=R h . Sin( ' l ' 4) 

mm 

Rc4= 12.23 


Rc5:=R b . Sin( ' 1,3) 

mm 

Rc5 = 14.122 

• 
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Computation of Plate Shape at Several Loads 


♦ 



Pj = 4.616x 10 3 N 
r := Rcl.. 110 


P 2 = 5.771 x 10 3 N 
r2:= Rc2.. 110 


P 3 = 7.694x 10 3 N 
r3:=Rc3.. 110 


4 

P 4 = 1.154x 10 N 
r4 := Rc4„ 110 


wl(r):= 


w2(r2) := 


w3(r3) := 


w4(r4) 


Pjln 

L r + 

a 

W - (si»( V l)) 4 l 


-L R + 

A 

^(R) 2 - (sin(xpl)) 4 _ 

j] 

[ 2 - n - N o 


mm 


— 

ho 

_El_ 

L r2 + 

./J2) 2 - (sin(vp2)) 4 J 

]11 

1r + 

_J(R) 2 - (sin(\|/2)) 4 ] 

I 


2' 7t 'N 0 


mm 


Pjln 


r3 + 

A 

(r3) 2 - (sin(\p3)) 4 _ 

]T1 

R + 

A 

f(R) 2 - (sin(\j/3)) 4 ] 

IJ 




2' 7t 'N 0 



mm 


P 4 ln 


_r4 + 

J(r4) 2 - (sin 

1-1 

]11 

R + 

A( R) 2 -(sin< 

W)"] 

IJ 



2.*-N 0 




mm 


P 5 = 1.539x 10 4 N 
r5:=Rc5.. 110 


w5(r5) := 


P 5 ln 


r5 + 

U 

(r5) 2 - (sin(v|/5)) 4 ]_ 


[r + 

[J 

^(R) 2 - (sinCvyS)) 4 _ 



2-n-N 0 


mm 


# 


87 




Displacement (mm) 


Plot of Analytic Solution for Plate Shape Up to Failure (Valid for r > $ 



r > r2,r3,r4,r5,r^ 
Plate Radius (mm) 


PUNCH SIZE: 


Rjj = 20mm 


r 0 > 


f 0 > 

4.616 


7.735 

5.771 


8.648 

7.694 

N 10 3 1^ = 

9.986 

11.541 


12.23 

15.388 


14.122 

^23.082> 


v. 17.296> 


P fail = 2.308 x 10 4 N 


Rc fail = 17.296mm 





Force (KN) 


Plot the Force vs. Displacement Relation: 







Displacement at Point C: 


Total Punch Displacement: 


8: 


° 1 

wl(Rcl) 

w2(Rc2) 

w3(Rc3) 

w4(Rc4) 

w5(Rc5) 


^ w fai] (R c 0 J 


5 plot := % 


(l -^(Vi)) 


-5; 


mm 



r 0 > 

9.499 

11.477 

14.635 

20.605 

26.311 

^37.637; 


6 max :=trunc 


V 


(1 - cos W) 


mm 


- w fail< Rcf ) 


Analytic Force vs. Displacement Curve 
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045" Thickness using 50 mm Radius Punch 


Entering Arguments for Analytic Solution 

Work Hardening Exponent and Coefficient: 

n := .22 C Q := 46010 6 -Pa 

Material Yield and Ultimate Strength Gy := 2 8600psi a u := 47600psi 


Plate Thickness 
Plate Clamped Radius 
Punch Radius 


t Q := ,045-in 
R := 110mm 


Rjj = 50 mm 


Compute Failure Parameters 


V c fail := (- 957 + -399n)-rad 


H 7 c fail ^ 


„ ... .33 .48 .52 

5 fail :=L41n R Rb 


Compute Stresses 


Flow Stress 


6 fail = 62.453mm 



+ °u) 


a 0 = 2.627 x 10 Pa 


^°° 7 > 


5 N 

N = 3.467x 10 - 
rn 


Generalized Membrane Stress 




Plate Predictions at Failure 


R := ) 10 


P fail :== . 2KN o R b( sin ( v ('c_fail)) 2 


P fai |= 8.l46x ]0 N 


R c_fail :=R b sin ('t' c _fail) 


R, 


Ref := 


c fail 


mm 


rr := Ref.. R 


& 


w 


fai 


K) := 


Pfail'N 


ff + 


R + 


J( r f) 2 - ( sin (^ c fail)) 4 


,j(R) 2 ~ ( sin (^c_fail)) 4 


2-rt-N 


o 


mm 


Establish Load Range in 5000 N Increments, Up to the Predicted Failure Load: 


Load: 

i := 1.. 7 


P := 


r n \ 


Mail 

5-N 

P fail 

4N 

P fail 

3N 

P fail 

2-N 


2P 


fail 


3N 
P fail 

V N J 


N 


P 0 :=0N 


Pi := 


r fail 


P 2 := 


P 3 := 


P 4 := 


P 5 := 


P fail 


P fail 


P fail 


2P fail 


P, = 1.629x 10 N 


P 2 = 2.037x 10 N 


P 3 = 2.715x 10 N 


P 4 = 4.073x 10 N 


P 5 = 5.431 x 10 N 
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Wrapping Angle 


( 


Vi = 


f 0 

22.753 
25.62 
29.954 
37.699 
44.92 
V 59.861 ) 


deg 


\j/0 := asin 


\|/2 := asin 


i|/4 := asin 


2.7fN 0 Rb 


2KN 0 Rb 


2-^N 0 -Rb 


\\)6 := asin 


[ r fail 

2-«'N 0 -Rb 


vj/j := asm 


2-"’N 0 -Rb 


vj/l := asin 


v[/3 := asin 


/ 

V> 


2-«-N 0 -Rb 


xp5 := asin 




2-nNoRb 


2' 7t 'N 0 -R b 




m 


Strain Components 

t 1 


e rr : = ln l 


V 


cos^i) 


e rr = 


e tt E rr 


' 0 > 


( 0 > 

0.081 


-0.081 

0.103 


-0.103 

0.143 

e tt = 

-0.143 

0.234 


-0.234 

0.345 


-0.345 

,0.689; 


v -0.689; 


Deformed Plate Thickness 

True Material Stress: 


t := t 0 cos(v|/j) 
°rr :=C o(~ In ( cos ( v i))) n 


m 


°rr = 


2.646x 10 


2.793 x 10 


3 x 10 


.8 


3.342 x 10 


3.64 x 10 


.8 


4.238 x 10 


Pa 
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To Find the Maximum Punch Displacement: 


R := 110mm 


Find the Amplitude: 



R 

— + 

_ R b 


( R ) 
UJ 

- (sin(v|/j)) 4 

sin(vj)(l + cos^j)) j 


Solve for Maximum Displacement using Geometric Relations 


< : = -%[l - cos( Vi ) + a-(sin(v|/j)) 2 J 




Find Radius of Contact Wrapping Point C: 


R c := R b -sin(x i / i ) 


*c = 


( 0 '' 
19.338 
21.62 
24.965 
30.576 
35.306 
V 43.241 j 


mm 
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Solve for Shape of plate, valid for < r < R 

R:= 110 


r 0 > 

16.292 

20.365 

P= 27.154 N-10 3 
40.731 
54.308 
l 81.462 J 


r 0 > 

22.753 

25.62 

vj/j = 29.954 deg 
37.699 
44.92 
i 59.861 ) 


Rc0:= R. s ' n 1'* ,0 l 
mm 

Rc0= 0 

Rcl^R b . ,i " ( * l) 

mm 

Rcl = 19.338 

Rc2:=R b . !ln( ' l ' 2) 

mm 

Rc2= 21.62 

Rc3 s R b . Sin(v3) 

mm 

Rc3 = 24.965 

Rc4:=R b Sin(,f4) 

mm 

Rc4 = 30.576 

Rc5 := Rl- 

mm 

Rc5 = 35.306 
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Computation of Plate Shape at Several Loads 


P, = 1.629x ]0 4 N 
r:= Rcl.. 110 


wl(r) := 


P, In 

Lr+Ll 

(r) 2 - (sin(v|/l)) 4 _ 


_[r + U 


]J 

W - (=in( v l))"] 


2^-N 0 



mm 


P 2 = 2.037x )0 4 N 
r2 := Rc2.. 110 


w2(r2) := 


— 

to 

.... 5“ 


_r2 + 

u 

(r2) 2 - (sin(v|/2)) 4 

]" 

R + 

J 


11 

^(R) 2 - (sin(v|/2)) 4 _ 




2-k.N 0 



mm 


P^ = 2.715 x 10 4 N 
r3:= Rc3.. 110 


w3(r3) := 



_r3 + 

( f 3) 2 - (sin(\|/3)) 4 ] 

P^-ln 

. L R + N 

f(R) 2 - (sin(\|/3)) 4 ]_ _ 



2.n-N 0 


mm 


P 4 = 4.073 x 10 4 N 
r4:=Rc4.. 110 


w4(r4) := 


p.-ln 

r4 + 

(r4) 2 - (sin(v|/4)) 4 ]_ 

i ^ in 

. _R + Li 

^(R) 2 - (sln(\^4)) 4 ] _ 



2-k-N 0 


mm 


P 5 = 5.431 x 10 4 N 
r5 :=Rc5.. 110 


P 5 ln 

[r5 + 

J (r5) 2 - (sin(\j/5)) 4 ]_ 


_Lr 4- 

J( R) 2 - (sin(iy5)) 4 __ 




2.tt.N 0 



w5(r5) := 




Displacement (mm) 


Plot of Analytic Solution for Plate Shape Up to Failure (Valid for r > g): 


Displacement vs. Plate Radius 



PUNCH SIZE: 


= 50mm 


< 0 > 


r 0 > 

16.292 


19.338 

20.365 


21.62 

27.154 

N-10 3 1^ = 

24.965 

40.731 


30.576 

54.308 


35.306 

^81.462; 


^43.24L 


P fail = 8 ,46x ,()4n 


Rc fail = 43.241mm 






Force (KN) 


Plot the Force vs. Displacement Relation: 


Displacement at Point C: Total Punch Displacement: 



f 0 ^ 


f 0 ^ 


wl(Rcl) 


16.893 


w2(Rc2) 

w3(Rc3) 

, .. ('-“W) . 

Splo, - Kb' mm -«i 6 pl«* 

20.125 

5j ;= 

25.165 


w4(Rc4) 

34.377 


w5(Rc5) 


42.928 


v w fail( Rcf ) y 


^59.814; 


8 max :=trunc 


%■ 


(l - cos(\|/6)) 


- w fail( Rcf ) 
























055" Thickness using 75 mm Radius Punch 


Entering Arguments for Analytic Solution 

Work Hardening Exponent and Coefficient: 

n := .22 C o :=46010 6 -Pa 

Material Yield and Ultimate Strength CTy : = 28600psi o u :=47600psi 


Plate Thickness 
Plate Clamped Radius 
Punch Radius 


t Q := ,055in 
R := 110 mm 


& 75 mm 


Compute Failure Parameters 

Vcjail := (- 957 + • 399n ) rad fail = 1045 

'2'2 40 C9 

5 fail : = 1 41 " R V Sfail = 77.111mm 


Compute Stresses 


Flow Stress 




) 


XT •- O 


a 0 = 2.627x 10 Pa 


Generalized Membrane Stress 





Plate Predictions at Failure 


R := 110 


P fail := _ 2 - 7tN o R b( sin (v c _fail)) 2 ] 


P fail= 1 493x 10 N 


*cjai\ :s Rb*sin( Vc _ fei |) 


R, 


Ref := 


c fail 


mm 


rf := Ref.. R 


9 


w fail( r f) : = 




_ r f + 

[J 

^f) 2 " ( sin ('t , c_fail)) 4 

IT 

tail 


R + 

0 

^R) 2 ~ ( sin ( v P c_fail)) 4 

]J 


2-7t-N„ 


mm 


Establish Load Range in 5000 N Increments, Up to the Predicted Failure Load: 


Load: 

i := 1.. 7 


P := 


f c\ \ 


r fail 

5-N 

P fail 

4-N 

P fail 

3N 

P fail 

2- N 

2P fail 

3- N 

P fail 

V ~n" j 


N 


P 0 :=0N 


P 1 : = 


p 2 := 


P fail 


P fail 


P 3 := 


P 4 := 


P 5 := 


"fail 


"fail 


2P fail 


Pj = 2.987x 10 N 


P 2 = 3.734x 10 N 


P 3 = 4.978x 10 N 


P 4 = 7.467x 10 N 


P 5 = 9.956x 10 N 
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Wrapping Angle 


( 




0 


\ 


22.753 
25.62 
29.954 
37.699 
44.92 
V 59.861 ) 


deg 


\[/j := asin 


( f P A 
2nN oRb 


vj/O := asin 


v[/2 := asin 


\|/4 := asin 


v|/6 := asin 


/ 

/ 

U 

f 

V 

/ 

V 


2-*-N 0 -Rb 


2-TtNoRb 


2-TtNoRb 


| Mail 
2-™-N 0 -Rb 


\|/1 := asin 


\j/3 := asin 


/ 

V 
/ 

V 
/ 


2nN 0 Rb 


2-w-NoRfc 


\j/5 := asin 


2-"-N 0 -Rb 





Strain Components 


( 1 




,n 1 

cos 

Vs) 

E tt ^ 

- £ rr 

V \ 

1 V ) 



' 0 > 


f 0 ' 

0.081 


-0.081 

0.103 


-0.103 

0.143 

8 tt = 

-0.143 

0.234 


-0.234 

0.345 


-0.345 

^ 0.689> 


^-0.689, 


Deformed Plate Thickness 

True Material Stress: 


t:= t 0 cos (vt/j) 

a rr :=C o(- |n ( cos ( v t / i))) n 


o 


rr 


0 1 

2.646x 10* 

2.793 x 10 8 

3 x ]0 8 

3.342 x 10 8 

3.64 x 10 8 

k 4.238x lO 8 ^ 


Pa 






m 


m 


m 
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To Find the Maximum Punch Displacement: 


R 110mm 


Find the Amplitude: 



*| ef 
_1 

-> 

( R ^ 

A 2 

- (sin(vi)) 4 

1 sin(vj)(l + cos( Vi )) I 


Solve for Maximum Displacement using Geometric Relations 


" cos (vj) + a-(sin(\|/j)) 2 ] 


Find Radius of Contact Wrapping Point C: 


R C := *b-si n (vi) 


R c = 


f 0 

29.007 
32.431 
37.448 
45.864 
52.959 
V 64.861,/ 
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Solve for Shape of plate, valid for < r < R 

R := 110 


P = 


r 0 > 


' ° 1 


r 0 N 

29.869 


0.397 


22.753 

37.337 


0.447 


25.62 

49.782 

N10 3 v|/j = 

0.523 

Vj = 

29.954 

74.673 


0.658 


37.699 

99.565 


0.784 


44.92 

J 49.347, 


J.045; 


^ 59.861 j 


deg 


Rc0:=R.- S,nW 

mm 

Rc0= 0 


Kcl»K b . ,in( ' l ' ,) 

mm 

Rcl = 29.007 

m 

Rc2 : .R h Sin(v2) 

mm 

Rc2= 32.431 


Rc3:=R b 5in( ' l ' 3) 

mm 

Rc3 = 37.448 

• 

Rc4*R b Sin( ' , ' 4) 

mm 

Rc4 = 45.864 


RcS : =R b . Sin( ' l ' 5) 

mm 

Rc5 = 52.959 

• 


m 


■# 


* 
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Computation of Plate Shape at Several Loads 




mm 



mm 



mm 



mm 
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Displacement (mm) 


Plot of Analytic Solution for Plate Shape Up to Failure (Valid for r > $ 



r,r2,r3,r4,r5,rf 
Plate Radius (mm) 


PUNCH SIZE: 


= 75mm 


f 



0 > 
29.869 
37.337 
49.782 
74.673 
99.565 


N-10 3 


V 149.347/ 




f 0 ^ 
29.007 
32.431 
37.448 
45.864 
52.959 


mm 


V 64.86lj 


P fail = 1 493x 10 5 N 


R, f a jj = 64.861mm 






Plot the Force vs. Displacement Relation: 


Displacement at Point C: Total Punch Displacement: 


f 0 > 


f o 

wl(Rcl) 


20.79 

w2(Rc2) 

(l -cos(yj)) 

24.502 

w3(Rc3) 

W-v' J -V «piot * 

30.165 

w4{Rc4) 


40.193 

w5(Rc5) 


49.228 

v w fai|C Rcf ), 


v 66.974; 


6 max :=trunc 


H- 


(l- 


cos(\|/6)) 

mm 


- w fail( Rcf ) 


Analytic Force vs. Displacement Curve 












Appendix G: Petalling Work Approximations 


Petalling Work Computation for .035” Sheet 
with 9.5mm Starter Hole and Four 3-mm Precracks 


Material Properties 

Gy = 23300psi g u = 45900psi 


a o := 


(°y + °u) 


Number of Petals 


N p : =4 


Crack Length 


aQ s 0 mm 


Petal Semi-Angle 

aj = 10mm a 5 = 50mm 

= 20mm a^ = 60mm 

a ^ = 30mm ay s 65 mm 

a 4 = 40mm ag s 70mm 

Moment Amplification Factor 

360 



0 := - deg 

2N p 

r| ^ 1 


Plate Thickness 

Petal Root Length 

t = .035 in 

b 0 := a o* sin ( 0 ) 



b| := a|-sin(o) 

b^ := a^-sin(o) 


b 2 := a 2' s * n ( 0 ) 

b 6 : = a 6 sin ( 0 ) 


b 3 := a 3 s i n (^) 

by := aysin(o) 


b 4 := a 4 sin(e) 

b 8 := a 8 -sin(e) 


Measured radius of Curvature 


p := 150mm 






Bending Moment 


2 V< b 0 

M o0 : =—7“ 

2a„t 2 -bi M . 2o o- ,2b 5 
o5= ^~“ 

2o| 2 b, 2c, o 1 b 6 

M o2 :=-^ M <*' = —T~ 


2o 0 -.-b 7 


2a t bo A/f . "o* 
w o j M„ 7 :=- 

M o3-=--- 07 4 


2a n t 2 b, M . 


Starter Hole Radius 


r Q := 4.75mm 


Cone Punch Horizontal Angle 

<|> := 30-deg 


Punch Vertical Speed 

V= ,a™ 

p min 

Petal Perpendicular Length 


:= aQ-cos(0) 

:= ajcos(o) 

l 5 := a^cos^) 

• 

:= a2-cos(o) 

1 6 := a 6 cos(e) 


:= a^-cos^) 

\-j .= ay COS ( 9 ) 


:= a 4 cos(e) 

*8 :=a 8 c °s( e ) 

• 


Rate of Petal Length change 

'dot : =tan(4.) v 
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Crack Tip Opening Distance (CTOD) 
CTO A 


CTOA = 10-deg 


5 c q := 2-(a 0 -sin(CTOA)) 

8 cl := 2-(a r sin(CTOA)) 

8 c5 2 ( 

5 c 2 := 2-(a2-sin(CTOA)) 

8 c6 :=2 ( 

5 c 3 := 2-(aj-sin(CTOA)) 

8 c7 := 2 ‘( 

8 c4 := 2-(a 4 -sin(CTOA)j 

8 c8 := 2 ( 


Punch Travel 

8 punch0 : =[ tan ( e ) ( l 0 + r o)] 
5 punch1 : =[‘an(e) (l, + rj] 
5 punch2 :=[tan(e) (l 2 + r Q )] 
8 punch3 :=[tan(e) (l 3 +r 0 )] 
8 punch4 : =[ tan ( 0 ) ( l 4 +r o)] 


5 sin(CTOA)) 

6 sin(CTOA)) 

8 sin(CTOA)) 

g-sin(CTOA)) 


8 punch5 : =[ tan ( 0 ) ( l 5+ r o)] 
8 punch6 • = [^ tan ( 0 )'('6 + r o)] 
8 punch7 - = j^tan{() (ly 4 r t) )j 

5 punch8 := [ tan ( 0 ) ('g + r 0 J] 





Rate of Work Done at Specified Increment 


W 


rateO 4 


W 


ratel = 


W rate 2 := 


W 


rate3 


M o0, 


"l0-tan(0)> 


(F) N \ 
°c0 

■~‘dof 

4-n- 


+ 3.84 


m 

- 

p ; 


< t j 

M ol , 


' lj tan(0)^ 


f 5 cl1 

‘dot* 

4rv 

-. 

+ 3.84| 

— 

m 

r 

\ 9 J 


1 t ) 

M o2, 


( l2’tan(0) N 


fy 

-’dot' 

4n- 

+ 3.84 

— 

m 

r 

V P ) 


1 t J 

M o3, 


r lytan(0)^ 


f y 

m dot* 

4-n- 


+ 3.84| 

— 

m 

_ 

l p J 


1 t J 


-4 


(sin(e)) 3 -(cos(e)) 1 



— (sin(e)) -(cos(e)) 


- ! 



- 4 
3 


- 1 


M 


W 


o4 


rate4 • 'dot 

m 


M 


W 


rate5 * 


o5 

m 




4*n* 


4’t] * 


+ 3.84 


ry 

t j 


M 


W 


06 


rate6 ■ 


m 


IVL 


**dot 


4 n 


W,..e7 = -f ! 'do.' 
m 


4-n 


l 4 -tan(e) ^ 


P J 


lrtan(o) 


( l 6» an ( 0 ) ^ 


v p ; 


Ijtan^)^ fS^V / n3 


(sin(e)) -(cos(e)) 


-4 

f | (sin(e)) 3 (cos(e)) 1 


— - 4 

'^ sin ^ e ^ 3 ( cos ( 0 )) 


(sin(©)) 3 (cos(e)) 



+ 3.84 


M 


W 


08 


rate8 •“ ‘dot 


m 


4-tj ■ 


V p J 


lolanCe) 


P J 


+ 3.84 


V t J 




(sin(e)) 3 (cos(e)) 1 


V t ) 


-4 

3 / /^w- 1 


“ I (sin(e)) *(cos(e)) 
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Total Membrane Work 


W 


mO' 


M 


(’o - r o) 


oO’ 


m 


3.84 


1 


[ t ] (f] 3 - (sin(9)) 3 


W ml : = 


W m2 : = 


w m3 : = 


M o r 


Ol-o) 


m 


1 




— t 


3.84 


P 


M 


fe-o 


o2* 


m 


3.84 


M 


o3‘~ 




m 


3.84 


V » J \X ) 

n3 

t j 

s \ 

°c3 

* J 


(sin(e)) 


^ Sin ^ 9 ^ 


— L 

3 


1 




•IfJ (-(e)) 


w 


m4* 


W m5 : = 


M 


('4 “ r o) 


o4' 


3.84 


J c4 

t ; 






— ] (sin(e)) 
t 


M 


o5‘~ 


Os-g 


m , 


1 


3.84 


n3 

^ t ) 


(f) <sin(e)) 


w m6~ 


M 


( ! 6- r o) 


06 ' 


m 


3.84 


W 


m7‘ 


M 0 r 


(*7 “ r o) 


r *_ c 6\ 3 ' " 3 

U J 

1 


— t 


f] (sin(9» 


3.84 


^) 3 (£L) 3 . W9 » 3 


W m8 : = 


M 


08 * 


Os-g 


m 


3.84 


/ x3 

V t j 


• (sin(G)) 

t 


— i 

3 


(cos(e)) 1 

(cos(e)) 1 

(cos(e)) 1 

(cos(e)) 1 

(cos(e)) 1 

(cos(e)) 1 

(cos(e)) 1 

(cos(e)) 1 

(cos(e)) 1 
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Total Bending Work 


' V b0-'VMo0('0- r o) 


2 tan(0) 


w bl - N P M ol ('l - r o) 2 ~ w bS » N p M o; .(l 5 - r 0 ) 2 

w b2 =■ vMb - 0 2 ^ W bs * n p . Mo6 .|, 6 - tfJS®. 

W b3 * N p Mo 3 ('3 - 'of ~ »b7 » Np Mor('7 ‘ 'of ~^ 


w b4*Np' M o4 0-|-'o) : 


2 pm 
2 tan(e) 
2 pm 


w b8- N pM o8 (l 8 -r 0 ) 


2 pm 
2 tan(Q) 
2 pm 


Membrane and Bending Work Totals 


W m o=ONm 

W b0 = 0N-m 


W m] = 0.641Nm 

W bJ = 4.788x 10" 

3 Nm 

W m 2= 6.54Nm 

W b2 = 1.568x 10" 

3 Nm 

W m3 = 13.648N m 

W b3 = 7.227x l(f 

3 Nm 

W m4 = 41.295Nm 

W b4 = 0.02N-m 


W m5 = 72.31 INm 

W b5 = 0.042N-m 


W m6 = 113.514Nm 

W b6 = 0.076Nm 


W m? = 141.606N m 

W b? = 0.098Nm 


W m8 = 165.581N m 

W b g = 0.125Nm 







Total Work 

W tO : = N p( W bO +W m o) 

W tl : = N p-( W bl + W ml) 

W t2 : = N p( W b2 +W m ^ 

W t3 : = N p-( W b3 +W m3) 

W t4 : =M W b4 +W m4) 
W t5 : =M W b5 + W m 5 ) 

W t6^= N p( W b6 +W m 6) 

W t7 : = N p( W b7 +W m7) 

W t8 ; = N p( W b8 +W m8) 


w t0 = oj 

W tl = 2.566J 
W t2 = 26.167J 
W t3 = 54.622J 
W t4 = 165.257J 
W t5 = 289.409J 
W t6 = 454.36J 
W t? = 566.816J 
W t8 = 662.824J 
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Estimation of Forces 


W 


F 0 := 2 ", 


F, := 2- 


tO 


^punchO 


W 


tl 


W 


F 5 :=2*> 


t5 


^unchl & punch5 


W 


F 2 :=2“ 


12 


W 


°punch2 
W* 


f 6 ^2: 


t6 


3 punch6 


F3 := 2*" 


t3 


W 


°punch3 


F 7 := 2 *“ 


t7 


°punch7 


W 


F 4 :=2- 


t4 


W 


F 8 := 2 " 


t8 


C 0 C 

°punch4 °punch8 


Force and Displacemnt Arrays 


V 


f 0 ^ 


f 0 ^ 




R 


434.109 

/ ~ \ 

F 1 


^punchl 




u 

F 2 


^punch2 


2.77 x 10 3 


11.821 





__ 3 


18.892 

F 3 


^punch3 


4.208x 10 

A 


25.963 

F 4 

6 := 

c 

°punch4 

P = 

1.001 X 10 

N 6 = 

33.034 

F 5 


c 

°punch5 


1.443x 10 4 


40.105 



r 


4 


47.176 



U punch6 


1.926x 10 









50.712 

F 7 


^punch7 


2.235 x 10 4 

! 

^54.247; 



^punch8 ) 


^ 2.444 x 10 4 y 



b 
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% 


Appendix H: LS-DYNA Input File 

$ .055" plate with 75 mm radius sphere 
$ 

^KEYWORD 

’TITLE 

Stamping 

TART 

part 

1 2 2 

$...+—1—-+—2-—+—3—+—4 —-+—5— + —6—+—7—+—8 
$ (2) CONTROL CARDS. 

$—+—i.—+ -—2-—+—3—+—4—+—5—+—^—+—7—+—8 
^define curve 
$ lcid (x-direction) 

1 

$ abscissa ordinate 
0.0 1000.0 
100.0 1000.0 
*control_termination 

$ ENDTIM ENDCYC DTMIN ENDNEG ENDMAS 
l.e-1 

*control_timestep 

$ DTINIT SCFT ISDO TSLIMT DTMS LCTM ERODE MS1ST 
.000 .900 0 

*CONTROL_SHELL 

$ WRPANG ITRIST IRNXX ISTUPD THEORY BWC MITER 
20.000 2 -1 1 2 2 1 

*CONTROL_HOURGLASS 
$ IHQ QH 

1 .100 

’ t CONTROL_BULK_VISCOSITY 
$ Q2 Q1 
1.500 .060 

*CONTROL_ENERGY 
$ HGEN RWEN SLNTEN RYLEN 

2 2 11 

TONTROLOUTPUT 

$ NPOPT NEECHO NREFUP IACCOP OPIFS IPNINT IKEDIT 
0 0 0 0 .000 0 100 

TATABASEEXTENTBINARY 

0 0 3 0 1 1 1 1 

0 0 0 0 0 0 

i—2—+--3—+--4 --+—5—+—6—+—7—+—8 
$ (4) DATABASE CONTROL CARDS FOR ASCII FILE 

$—|-—1-1- —2-1 —3—+—4—+—5—+—6—+—7——8 


% 

% 


* 
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"DATABASERWFORC 
l.e -4 

*DATABASE_GLSTAT 

l.e -4 

*DATABASE_MATSUM 

l.e -4 

$—i-—i f —2 1 —3 1 —4 v — 5 —-i—6—(- — 7 —i—8 

$ ( 5 ) DATABASE CONTROL CARDS FOR BINARY FILE 

$--+--1 --+-— 2 —+— 3—+—4 --+— 5 —+--6—+— 7 —+—8 
*DATABASE_BINARY_D 3 PLOT 
$ DT/CYCL LCDT NOBEAM 
l.e -3 

$-+— 1 —+— 2-—+—3 --+— 4 —+— 5 —+-6--+— 7 —+—8 
$ RIGID_SPHERE 

— i — + — 2-—+ — 3 — + — 4 —+— 5 —+—6—+ — 7 — + —8 
*rigidwall_geornetric_spheremotion 
$ NED/NSID 
0 

S xt yt zt xh yh zh fric 

0 . 100 . 75.7 0 . 100 . 0 . 0.05 

$ radius 
75.0 

$ 

$ LCID OPT VX VY VZ 
1 0 0. 0. -1.0 

$—i-—i f-— 2 —+— 3—+—4 h— 5 —+—6—h 7 — + —8 

$ BOUNDARY CONDITION 


$-1-1-1- 

— 2—+—3 -h- 

-4 -|- 5 __ 

--+—6— 

-+— 7 —+■ 

—8 


^BOUNDARY SPC SET 






Sclamped 







$ NID/NSID 

CED DOFX 

DOFY 

DOFZ 

DOFRX 

DOFRY 

DOFRZ 

1 

1 1 1 

1 1 

1 




^BOUNDARY SPC SET 






$symy 







$ NID/NSID 

CED DOFX 

DOFY 

DOFZ 

DOFRX 

DOFRY 

DOFRZ 

2 

1 

1 

1 




^'BOUNDARY SPC SET 






$symx 







$ NID/NSID 

CID DOFX 

DOFY 

DOFZ 

DOFRX 

DOFRY 

DOFRZ 

3 

1 

1 

1 




"BOUNDARY SPC SET 






Scenter 







$ NID/NSID 

CID DOFX 

DOFY 

DOFZ 

DOFRX 

DOFRY 

DOFRZ 

4 

1 1 

1 1 

1 




$—i-1-1-. 

— 2 —+— 3 —— 

-- 4 —+— 5 - 

-+—6— 

•+— 7 — + - 

-8 



=: MAT_PIECEWISE_LINEAR_PLASTICITY 
$ MID RO E PR SIGY ETAN 
2 7 . 800 E -09 210 .E +03 0.33 


EPPF TDEL 




S C P LCSS LCSR 

$ EPS1 EPS2 EPS3 EPS4 EPS5 EPS6 EPS7 EPS8 

0 .02 .04 .06 .1 .5 .75 1 

$ ESI ES2 ES3 ES4 ES5 ES6 ES7 ES8 

240 260 295 330 695 920 1160 2000 

*SECTION_SHELL 
SHMNAME PROPS lSectShlll 

$ SID ELFORM SHRF NIP PROPT QR/IRID ICOMP SETYP 
2 2 0.833 5 

$ T1 T2 T3 T4 NLOC 

1.397 1.397 1.397 1.397 

$—+—i—+—2-—+—3—+—4—+—5—+—6—+—7—+—8 
S GEOMETRY 

$—h—1 1—2—i—3 —i 4 —i—5 —i—6—i 7 i 8 

*NODE 
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Appendix I: Strain Field Worksheet 


Strain Field Computation Worksheet 

ROW 0 


u i i_i_ptoo = 0 u n_ 2 _ptoo = 0 


^11 _pt00 - 1 0 § y 22_pt00 = 1 0 

u 22_l_pt00 = 0 u 22_2_pt00 = 0 


^MljrtOO := U 1 l_2_pt00 " U 1 l_l_pt00 

8u 22_pt00 : - u 22_2_pt00 

^lljtOO 

5u 22_pt00 

E ll_pt00~ , 

“1 l_pt00 

£ 22_pt00 - x 

5 y22_pt00 

lIT^H-PtOo) f 5u 22_pt00 Yl 



' 12_pt00 • 


2 Ll 5y 22jpt00 J 


V 5x ] 1 _ptOO JJ 


=-ll_ptOO 

E 22_pt00 
E 12_pt00 


= 0 
= 0 
= 0 


Next Point 


• 

u ll_l_ptl0= 8 
^1 l_ptl0 = - 8 
u 22_l_pt10= 0 

u ll_2_ptl0 = 10 
^22_ptl0= 8 
u 22_2_ptl0- 0 


• 

Su ll_ptl0 : - u ll_2_ptl0~ u ll_l_pt10 

® u 22__ptl0 u 22_2_j)tl0“ u 22_l_ptl0 


^lljtlO 
e n_ptio - , 

6x ll_ptlO 

6u 22_ptlO 
£ 22 ptlO — , 

S >22ptl0 

• 

1 

e 12_ptl 0" 2 

<5u ll_ptlO^| i ^ u 22_ptlo'' 
v 5 y22_ptloJ V 6x ]l_ptl0; 

e l l_ptl0 = 0 25 
E 22_ptl0 = 0 
E ]2_ptl0 = 0125 
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Epsilon -n Strain Field: 


£ 


11 := 


E 1 l_pt06 E 1 l_ptl6 e 11_pt26 £ 1 l_pt36 

e ll_pt05 e l l_ptl5 e 1 l_pt25 E 1 l_pt35 
e ll_pt04 e ll_ptl4 e ll_pt24 e ll_pt34 
E 1 l_pt03 e l l_pt13 E 1 l_pt23 E 1 l_pt33 
e l l_pt02 e 11 __pt12 e l l_pt22 e 11 _pt32 
E ll_pt01 E ll_ptll E 1 l_pt21 £ 1 l_pt31 
< E ll_ptOO E ll_ptlO e l l_pt20 E 1 l_pt30 


E 1 l_pt46 e 1 l_pt56 E 1 l_pt66 E 11 jjt76^ 
e l l_pt45 e ll_pt55 E lljpt65 e ll_pt75 
e l l_pt44 E 1 l_pt54 E ll_pt64 e ll_pt74 
E 1 l_pt43 E 1 l_pt53 e ll_pt63 E 1 l_pt73 

E 1 l_pt42 E 1 l__pt52 E ll_pt62 e l 1_pt72 

£ 1 l_pt41 E ll_pt51 e l l_pt61 E 11_pt71 

E llj3t40 e ll_pt50 e ll_pt60 E 1 l_pt70> 



f 0 

0.1 

0.15 

0 

0 

0.2 

0.4 

0.25 ^ 


0.5 

0.1 

0.2 

0.05 

0.2 

0.3 

0.4 

0.3 


-0.222 0.111 

0.2 

0.1 

0.2 

0.4 

0.4 

0.5 

E 11 - 

-0.25 

0.222 

0.222 0.278 0.333 

0.444 

0.722 0.579 


0 

0.111 

0.167 0.222 0.333 

0.444 

1.056 0.842 


0 

0.111 

0.25 

0.25 

0.5 

0.5 

0.625 

0.941 


< o 

0.25 

0.25 

0.375 0.571 

0.75 

0.5 

0.75 , 



Epsilon ?? Strain Field: 


£ 22 : " 


s 22_pt06 

£ 22_pt05 

£ 22_pt04 

s 22_pt03 

s 22_pt02 

£ 22_pt01 

V £ 22_pt00 


£ 22_ptl6 
e 22_ptl 5 
£ 22_ptl4 
e 22_ptl3 
£ 22_ptl2 

£ 22_ptl 1 
£ 


e 22_pt26 

e 22_pt25 

E 22_pt24 

e 22_pt23 

£ 22_pt22 

£ 22_pt21 
E' 


22_ptl0 t 22_pt20 


£ 22_pt36 
£ 22_pt35 
E 22_pt34 
£ 22_pt33 
e 22_pt32 
£ 22_pt31 
£ 22_pt30 


e 22_pt46 

£ 22_pt45 

£ 22_pt44 

e 22_pt43 

e 22_pt42 

E 22_pt41 


t 22_pt40 t 22_pt50 


E 22_pt56 

e 22_pt55 

e 22_pt54 

e 22_pt53 

£ 22_pt52 

e 22_pt51 
£' 


e 22_pt66 

E 22_pt65 

e 22_pt64 

e 22_pt63 

E 22_pt62 

E 22_pt61 

E 22_pt60 


e 22_pt76 >1 

e 22_pt75 

£ 22_pt74 

e 22_pt73 

£ 22_pt72 

e 22_pt71 

£ 22_pt70 J 


z 22~ 


0.4 

0.5 

0.55 

0.6 

0.7 

0.7 

0.8 

0.8 

0 

0.5 

0.55 

0.65 

0.7 

0.9 

1 

3 

0.2 

0.3 

0.45 

0.45 

1.7 

1 

1 

1.05 

0.2 

0.111 

0.222 

0.333 

0.389 

0.944 

1 

1 

0 

0 

0.167 

0.222 

0.222 

0.222 

3.1 

3.15 

0 

0 

0 

0 

0.1 

0.15 

0.2 

0.25 

0 

0 

0 

0 

0 

0.125 

0.187 

0.25 
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Epsilon 1? Strain Field: 


£ 12 : ~ 


<e 12_pt06 £ 12_ptl6 £ 12_pt26 £ 12_pt36 e 12_pt46 e 12__pt56 
£ 12_pt05 £ 12_pt15 £ 12_pt25 £ 12_pt35 £ 12_pt45 e 12_pt55 
c 12_pt04 £ 12_ptl4 e 12_pt24 e 12_pt34 e 12_pt44 e 12_pt54 
£ 12_pt03 E 12 pt 13 e 12_pt23 £ 12_pt33 £ 12_pt43 e 12_pt53 
£ 12_pt02 e 12_ptl2 e 12_pt22 £ 12_pt32 £ 12_pt42 £ 12_pt52 
£ 12_pt01 £ 12_ptll e 12_pt21 e 12_pt31 e 12_pt4] e 12_pt51 
^ e 12_pt00 e 12_pt 10 e 12_pt20 e 12_pt30 £ 12jpt40 £ 12_pt50 


£ 12_pt66 

£ 12_pt65 
E 12_pt64 
e 12_pt63 

e 12_pt62 

e 12_pt61 

e 12_pt60 


f 0.2 

0.3 

0.35 

0.3 

0.35 

0.45 

0.6 

0.525^ 

0.25 

0.3 

0.375 

0.35 

0.45 

0.6 

0.7 

0.65 

0.011 

0.217 

0.325 

0.275 

0.95 

0.7 

0.7 

0.775 

0.025 

0.167 

0.222 

0.306 

0.361 

0.694 

0.881 

0.801 

0 

0.063 

0.167 

0.222 

0.278 

0.333 

1.086 

1.005 

0 

0.063 

0.125 

0.125 

0.263 

0.294 

0.375 

0.547 

< o 

0.125 

0.125 

0.187 

0.25 

0.437 

0.344 

0.5 > 


Rotating the strain field in order to align coordinate sysytem with crack: 

e 22 

£r m £ 22' cos < 22 de 8) £ ^ := — (22 deg) £ 4n := £ 12’C°s(22 deg) 


£ 12_pt76 > 

E 12_pt75 

£ 12_pt74 

e 12_pt73 

e 12_pt72 

e 12_pt71 

e 12_pt70 > 




Final Crack-Oriented Strain Fields: 


-Tin 


' 0.371 0.464 0.51 0.556 0.649 0.649 0.742 0.742^ 
0 0.464 0.51 0.603 0.649 0.834 0.927 0.927 

0.185 0.278 0.417 0.417 1.576 0.927 0.927 0.974 
0.185 0.103 0.206 0.309 0.361 0.876 0.927 0.927 
0 0 0.155 0.206 0.206 0.206 1.02 1.066 

0 0 0 0 0.093 0.139 0.185 0.232 

< 0 0 0 0 0 0.116 0.174 0.232 \j 


E tt = 


( 0.431 0.539 0.593 0.647 0.755 0.755 0.863 0.863^ 

0 0.539 0.593 0.701 0.755 0.971 1.079 1.079 

0.216 0.324 0.485 0.485 1.834 1.079 1.079 1.132 

0.216 0.12 0.24 0.36 0.419 1.019 1.079 1.079 

0 0 0.18 0.24 0.24 0.24 1.186 1.24 

0 0 0 0 0.108 0.162 0.216 0.27 

K. 0 0 0 0 0 0.135 0.202 0.27, 



e 4n = 


f 0.185 0.278 0.325 0.278 0.325 0.417 0.556 0.487^ 

0.232 0.278 0.348 0.325 0.417 0.556 0.649 0.603 

0.01 0.201 0.301 0.255 0.881 0.649 0.649 0.719 

0.023 0.155 0.206 0.283 0.335 0.644 0.816 0.743 

0 0.058 0.155 0.206 0.258 0.309 1.007 0.932 

0 0.058 0.116 0.116 0.243 0.272 0.348 0.507 

< 0 0.116 0.116 0.174 0.232 0.406 0.319 0.464, 
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